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Abstract 

The  dual  basic  tasks  of  evaluating  ship  waves  at  the  free  surface  and  of  removing  unwanted  short  waves  are 
considered  within  the  framework  of  the  ‘free-surface  Green  function  potential  flow  theory’,  based  on  a  Green 
function  that  satisfies  the  radiation  condition  and  the  Kelvin-Miclrell  linearized  boundary  condition  at  the  free 
surface.  A  practical  approach  based  on  parabolic  extrapolation  within  an  extrapolation  layer  bordering  the 
free  surface  is  used.  The  height  of  the  extrapolation  layer  is  defined  explicitly  via  simple  analytical  relations 
in  terms  of  the  Froude  number  and  the  slenderness  of  the  ship  hull,  and  varies  from  the  bow  to  the  stern.  The 
bow-to-stern  variation  is  an  important  ingredient  that  accounts  for  the  fact  that  waves  along  the  ship  hull  aft 
of  the  bow  wave  differ  from  the  bow  wave.  Indeed,  a  ship  bow  wave  is  significantly  higher  and  shorter  than 
waves  aft  of  the  bow  wave,  is  affected  by  nearfield  effects  related  to  the  rapid  variation  of  the  hull  geometry 
at  a  ship  bow,  and  consequently  contains  more  short  wave  components. 


1  Introduction 


An  offshore  structure  (or  any  other  floating  body  without  mean  forward  speed)  in  ambient  time-harmonic 
(regular)  waves  with  frequency  u>  generates  a  system  of  (diffracted-radiated)  waves  that  all  have  the  same 
frequency  ui  and  wavelength  A  =  2irg/uj2.  Such  a  system  of  waves,  with  a  single  discrete  wavelength,  differs 
in  a  fundamental  respect  from  the  waves  —  with  a  continuous  spectrum  of  wavelengths  —  that  are  generated 
by  a  ship  advancing  in  calm  water  or  in  ambient  regular  waves.  In  particular,  for  the  simplest  case  of  a  ship 
advancing  with  forward  speed  Vs  in  calm  water  considered  here,  the  ship  creates  waves  with  a  spectrum  of 
wavelengths  A  within  the  range  0  <  A  <  2tt V2/g.  Thus,  ship  waves  defined  within  the  classical  framework 
of  potential  flow  theory  include  very  short  waves  that  can  be  significantly  affected  by  surface  tension  and 
viscosity,  and  consequently  are  physically  unrealistic.  The  spectrum  of  ship  waves  may  also  include  waves  that 
are  not  appreciably  affected  by  surface  tension  and  viscosity  but  are  short  with  respect  to  the  ship  length  Ls , 
and  consequently  may  have  limited  effect  on  the  ship  drag,  sinkage  and  trim.  There  is  then  a  practical  need 
for  eliminating  short  waves  in  the  spectrum  of  waves  generated  by  ships.  Indeed,  an  effective  and  practical 
method  for  filtering  short  gravity  waves  is  an  important  ingredient  of  any  numerical  method  for  computing  ship 
waves.  This  basic  issue  is  considered  here  within  the  framework  of  the  ‘free-surface  Green  function  potential 
flow  theory’,  based  on  a  Green  function  that  satisfies  the  radiation  condition  and  the  Kelvin- Michell  linearized 
boundary  condition  at  the  free  surface. 

Thus,  we  consider  linear  potential  flow  about  a  ship  hull  of  length  Ls  that  steadily  advances  at  speed  Vs 
along  a  straight  path  in  calm  water  of  effectively  infinite  depth  and  lateral  extent.  The  flow  about  the  ship 
hull  is  observed  from  a  righthanded  moving  system  of  orthogonal  coordinates  X  =  (X,Y,Z)  attached  to  the 
ship,  and  thus  appears  steady  with  flow  velocity  given  by  the  sum  of  an  apparent  uniform  current  (—14, 0, 0) 
opposing  the  ship  speed  Vs  and  the  (disturbance)  flow  velocity  U  =  (U,V,W)  due  to  the  ship.  The  X  axis  is 
chosen  along  the  path  of  the  ship  and  points  toward  the  ship  bow.  The  Z  axis  is  vertical  and  points  upward, 
with  the  mean  (undisturbed)  free  surface  taken  as  the  plane  Z  =  0 .  The  length  Ls  and  the  speed  Vs  of  the 
ship  are  used  to  define  nondimensional  coordinates  x  =  X/Ls  and  flow  potential  <f>  =  &/(VsLs). 

We  define  the  usual  Froude  number  F=VS  f\JgLs  where  g  is  the  acceleration  of  gravity.  We  also  define 
the  Froude  number  Fbd  =  U  !\J gLBD  based  on  a  transverse  dimension  LBD  of  the  ship  hull  that  is  chosen 
as  Lbd  =  BD/(B/2  +  D ).  We  have  LBD  <  B  and  LBD  <  2D.  The  Foude  numbers  F  and  Fbd  based  on  the 
length  Ls  or  the  transverse  dimension  LBD  of  the  ship  are  related  as 


Fbd  = 


where  a 


(1) 
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As  already  noted,  the  flow  about  the  ship  is  considered  within  the  framework  of  the  ‘free-surface  Green 
function  method’  based  on  a  Green  function  G(x;x)  that  satisfies  the  radiation  condition  and  the  Kelvin- 
Michcll  linearized  boundary  condition  at  the  free  surface.  The  points  x  =  ( x,y,z )  and  x  =  ( x,y,z )  stand 
for  ‘boundary  points’  located  on  the  ship  hull  surface  T,H ,  and  ‘flow-field  points’  located  on  or  in  the  flow 
region  outside  Y,H . 

Within  this  approach,  the  flow  potential  </>  =  c/(x)  at  a  flow  field  point  x  can  be  expressed  as  <fi  =  (j)L+  4>w 
where  <f>L  represents  a  nonoscillatory  local  flow  component  and  <pw  defines  the  waves  generated  by  the  ship. 
These  local  and  wave  components  are  associated  with  the  decomposition  G  =  L  +  W  of  the  Green  function 
G  into  a  local  flow  component  L  and  a  wave  component  W,  as  in  e.g.  [1] .  The  local  flow  potential  <j>L  is  not 
considered  here. 

Thus,  we  only  consider  the  wave  potential  (j>w .  Within  the  free-surface  Green  function  approach,  the 
wave  potential  (f>w  =  0M/(x)  at  a  flow  field  point  x  =  (x,y,z<  0)  is  expressed  as  a  Fourier  superposition  of 
elementary  plane  waves  E .  Specifically,  <f>w  is  given  by 

i  n  oo  _  _  _ 

<j>w=-Sm  dkSE  with  E  =  eV+k2)z/F2  +  iVi+W  (x  +  ky)/F2  (2) 

*  J- oo 

e.g.  [2,3] .  The  amplitude  S  =  S(k;x)  of  the  elementary  waves  E  in  the  Fourier  representation  (2),  called 
wave-spectrum  function  or  Kochin  function,  is  also  given  by  a  Fourier  superposition  of  elementary  plane  waves. 
Specifically,  the  wave-spectrum  function  S  in  (2)  is  given  by  a  distribution  of  elementary  waves  E  over  the 
portion  T,H  of  the  mean  wetted  ship  hull  surface  T,H  that  is  defined  by  x  <  x.  We  then  have 

S  =  — ^  f  daAE  with  E  =  e(i+k2)z/F2 -iVi+V  (x+ky)/F2  (3) 

F  Jy.h 

Here,  da  =  da(x)  is  the  differential  element  of  area  at  a  point  x  of  the  ship  hull  surface  T,H . 

An  important  particular  case  of  the  generic  Fourier-Kochin  representation  (2)-(3)  of  ship  waves  is  the  wave 
potential  </>^  associated  with  the  Hogner  slender-ship  potential  4>h  given  in  [4] .  In  this  special  case,  the 
amplitude  A  of  the  elementary  wave  E  in  the  wave-spectrum  function  S  is  given  by 

A  =  Ah  with  AH  =  nx  (4) 

Here,  n  =  (nx,nv,nz)  is  a  unit  vector  that  is  normal  to  T,H  at  x  and  points  outside  T,H ,  i.e.  into  the  water. 
The  Hogner  slender-ship  approximation  (4),  defined  explicitly  in  terms  of  the  hull  geometry,  is  useful  for  many 
practical  applications.  Furthermore,  this  explicit  flow  approximation  is  a  major  element  of  the  Neumann- 
Michcll  (NM)  theory  of  ship  waves  given  in  [3]  and  considered  further  on.  Indeed,  the  NM  theory  provides  a 
correction  of  the  Hogner  slender-ship  potential  (f>H  given  in  [4] . 

The  wavelength  A  of  the  elementary  waves  E  and  E  in  (2)  and  (3)  is  given  by 

0  <  A  =  2ttF7(1  +  k2)  <  2nF2  =  X0  (5a) 

where  Ao  is  the  wavelength  of  the  transverse  waves  generated  by  a  ship  along  its  track.  Transverse  and 
divergent  ship  waves  correspond  to  \k\  <  l/\/2  «  0.71  and  1  /x/2  <  |Aj,  respectively;  e.g.  [5].  The  ratio 

A/Ao  =  1/(1  +  k2)  (5b) 

is  equal  to  1/2  for  k  =  1,  1/3  for  k  =  1/5  for  k  =  2,  1/10  for  k  =  3,  and  approximately  6%  or  4%  for 

k  =  4  or  k  =  5.  Thus,  values  of  |fc|  greater  than  3  correspond  to  waves  that  are  significantly  shorter  than  the 
longest  waves  created  by  a  ship. 

The  continuous  spectrum  of  waves  generated  by  a  ship  can  be  usefully  divided  into  ‘long  waves’  associated 
with  the  range  —kiong  <  k  <  kiong  where  the  cutoff  wavenumber  kiong  can  reasonably  be  taken  as  3  or  even  2, 
and  ‘short  waves’  that  correspond  to  kiong  <  |/c|.  These  short  waves  can  be  further  divided  into  ‘short  gravity 
waves’  that  are  too  long  to  be  significantly  affected  by  surface  tension  and  viscous  effects,  and  ‘very  short 
waves’  for  which  surface  tension  and  viscous  effects  cannot  be  ignored.  An  elegant  physics-based  theory  that 
accounts  for  the  influence  of  surface  tension  and  viscosity  on  gravity  waves  generated  by  a  ship  hull  advancing 
in  calm  water  or  in  ambient  time-harmonic  waves  is  expounded  in  [6-9] .  Surface  tension  and  viscosity  are 
ignored  here  for  simplicity,  and  because  the  cutoff  wavenumber  kiong  is  expected  to  correspond  to  relatively 


long  waves  not  significantly  affected  by  surface  tension  and  viscosity  (although  that  may  not  always  be  the 
case  at  model  scale). 

For  a  fully  submerged  body,  for  which  2  <  —S  where  0  <  5  is  the  distance  between  the  free  surface  plane 
z  =  0  and  the  highest  point  of  the  submerged  body  surface  (point  nearest  the  mean  free  surface),  we  have 
\E\  <  e_(-1  +  fc  )S/F  .  Thus,  the  elementary  wave  function  E  and  consequently  the  spectrum  function  S  decay 
exponentially  as  k  — >  ±00  for  a  fully  submerged  body.  The  Fourier  integral  (2)  therefore  converges  for  every 
value  of  z  <  0,  and  does  not  contain  short  waves,  in  this  special  case. 

However,  the  situation  is  different,  and  more  complicated,  for  a  surface-piercing  ship  hull  T,H ,  for  which 
z  <  0  and  the  function  E  does  not  decay  exponentially  as  k  — »  ±00.  If  z  <  —h^<  0,  i.e.  for  flow  field  points  x 
at  some  distance  below  the  mean  free  surface  z  =  0,  the  exponential  function  E  in  (2)  decays  exponentially  in 
the  limit  k  — >  ±00  and  the  Fourier  integral  (2)  can  be  evaluated  accurately.  In  practice,  the  infinite  limits  of 
integration  in  (2)  are  replaced  by  finite  limits  ±k0 a.  The  function  E  is  smaller  than  0.7%  if  {l  +  k2)z/F2  <  —  5, 
and  (2)  can  be  evaluated  accurately  for  z  <  —hoo  with 

h00  =  bF2/{l  +  kl0)  (6) 

Convergence  of  the  Fourier  integral  (2)  is  not  a-priori  obvious  for  —hoc  <  ~z  <  0,  i.e.  for  flow  field  points  x  in 
the  vicinity  of  the  free  surface  z  =  0.  Indeed,  we  have  E  =  1  if  z  =  0 ,  for  every  value  of  the  Fourier  variable  k . 

Furthermore,  the  relation  (6)  yields  k^  — >  00  as  h0 0  — »  0,  and  therefore  implies  that  the  wave  potential 
<j)W  includes  very  short  waves  in  this  limit.  However,  very  short  gravity  waves  are  unrealistic  because  surface 
tension  and  viscous  effects,  ignored  in  (2)  and  (3),  cannot  be  neglected  in  the  short-wave  limit  k  — >  00,  and 
because  short  waves  that  correspond  to  kiong  <  |  fc  are  of  limited  interest  for  most  practical  applications.  Thus, 
robust  evaluation  of  the  wave  integral  (2)  for  —hoo  <  z  <  0  is  a  nontrivial  basic  issue  for  the  computation  of 
ship  waves  within  the  framework  of  the  free-surface  Green  function  approach. 

A  practical  way  of  evaluating  the  Fourier  integral  (2)  at  and  near  the  mean  free  surface  z  =  0  is  then 
required.  To  this  end,  the  Fourier  representation  (2)  can  be  modified  as 

j>w=  -Qmf  °°dkKSE  (7) 

n  J-koo 

Selection  of  an  appropriate  finite  limit  of  integration  (cut-off  wavenumber)  koo ,  and  an  effective  short-wave 
filter  function  A  are  important  elements  of  the  free-surface  Green  function  theory  of  ship  waves  as  already 
noted. 


2  Parabolic  extrapolation 

A  simple  way  of  evaluating  the  Fourier  integral  (2)  for  flow-field  points  x  within  a  layer  —h  <  z  <  0  near 
the  mean  free-surface  plane  z  =  0  is  via  extrapolation  of  the  wave  potential  (f>w  computed  at  points  x  in  the 
region  z  <  —h.  In  particular,  the  wave  potential  <f>^  at  a  point  Xo  of  the  mean  free  surface,  i.e.  with  zq  =  0 , 
can  be  determined  via  parabolic  extrapolation  of  the  wave  potentials  <j>Y ,  <t>2  ,  fiY  computed  at  three  points 
xi ,  X2  ,  X3  with  Z3  <  Z2  <  Zi  <  —h .  Specifically,  the  parabolic  approximation 

~w  _  (z  -  z2)(z  -  Z^)4>Y  +  (z  -Z!){z  -Z3)4>f  +  {z  -  Zi)(z  -  Z2)4>3 
(zi  -  22)(?'l  -  Z3)  (z2  -  Zi){Z2  -  z3)  (Z3-Zl)(z3-Z2) 

with  z  =  0  yields  the  extrapolation  formula 

<!>q  ~  3 ( —  </>2V  )  +  (f> w  if  zi  =  —  h ,  z2  =  —2 h,z3  =  —3 h 

This  parabolic  extrapolation  formula  amounts  to  replacing  the  exponential  function  E  =  e  ( 1+k  ^  Z!F  in  (2) 
by  the  function 

E0  =  Eh(3-3Eh  +  E2)  with  Eh  =  e^1+k2^/F2 

for  z  =  0.  We  have  Eh  — >■  0  and  Eq  ~  3  Eh  as  k  — >  ±00  .  Thus,  while  E  =  1  for  z  =  0 ,  Eq  vanishes  in  the 
limit  k  — >  ±00  .  The  function  .E*  defined  as 

E*  =  E0  (1+  z/h)  —  Ehz/h  =  Eh  [3  —  3E^  +  E\  +  (2  —  3Eh  +  El)z/h]  with  —  h  <  z  <  0 
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Figure  1:  Variations  of  the  midship  and  bow- wave  extrapolation  heights  hm  and  hb  (left  side)  and  of  the 
cutoff  wavenumbers  koo  and  kiaT  (right  side)  within  the  range  0  <  F  <  0.8  for  Ft  =  0.4,  aH  =  0.06,  £b  =  0.04, 
£m  =  0.06  and  C™  =  0.7. 


can  be  substituted  for  the  function  E  within  the  layer  —h<z<  0  near  the  free  surface  z  =  0.  The  wave 
potential  ipw  is  not  affected  for  z  <  —h  if  the  exponential  function  E  in  the  Fourier  integral  (2)  is  replaced 
by  the  function  E*  =  AE  where  the  function  A  is  defined  as 

!A  =  1  for  z  <  —h  'j 

A  =  e-  d  +  fc2)  (^+^)/^2  [3  —  3Eh  +  El  +  (2- 3Eh  +  E2)z/h]  for-h<~<0>  (8) 

with  Eh  =  e~(l  +  k2^/F2  ) 


The  replacement  of  the  function  E  in  (2)  by  the  function  E*  =  AE  amounts  to  multiplying  the  integrand 
of  the  Fourier  integral  (2)  by  the  function  A,  as  in  (7).  This  function  provides  a  reasonable  extension  of  the 
wave  potential  (j>w  (which  is  unaffected  for  z  <  —  h)  up  to  the  free  surface  z  =  0,  and  ensures  convergence  of 
the  Fourier  integral  (7).  The  function  A  and  the  finite  limits  of  integration  ±fcoo  act  as  a  filter  that  eliminates 
short  waves  in  (2). 

The  short-wave  filter  (8)  depends  on  z/F2  and  h/F2.  Furthermore,  the  ‘extrapolation  height’  h  of  the 
‘extrapolation  layer’  —h  <  z  <  0  can  (indeed,  should)  be  chosen  to  vary  with  x,  as  shown  below.  The  filter 
(8)  is  more  elaborate  than  the  simpler  filter  used  in  [2,10],  which  is  independent  of  x  and  z.  The  filter  (8)  is 
also  more  effective,  notably  because  an  extrapolation  height  h  that  varies  along  the  length  of  a  ship  hull  is 
useful  as  now  shown. 


3  Basic  extrapolation  heights  and  related  cutoff  wavenumbers 

Waves  along  a  ship  hull  at  midship  —  at  some  distance  away  from  the  ship  bow  and  stern  where  the  hull 
geometry  varies  rapidly  —  are  largely  determined  by  the  ship  speed.  Indeed,  ‘midship  waves’  have  a  wavelength 
approximately  equal  to  the  (nondimensional)  wavelength  Ao  =  2ttF2  of  the  longest  waves  generated  by  a 
ship  along  its  track.  The  fundamental  wavelength  Ao  is  approximately  equal  to  1  (the  nondimensional  ship 
length)  for  F  ss  0.4.  The  extrapolation  height  h  for  the  waves  along  a  ship  hull  at  midship  (aft  of  the 
bow  wave  and  ahead  of  the  ship  stern)  can  then  be  taken  as  a  small  fraction  0  <  £m  <C  1  of  the  longest 
wavelength  Ao,  i.e.  hm  =  £m2nF2.  This  extrapolation  height  increases  rapidly  as  the  Froude  number  F 
increases,  and  ultimately  becomes  meaningless.  Indeed,  at  high  Froude  numbers,  the  waves  generated  by  a 
ship  are  increasingly  influenced  by  the  variation  of  the  hull  geometry,  which  is  primarily  determined  by  the 
hull  slenderness  aH  defined  by  (1).  The  extrapolation  height  hm  can  then  be  assumed  to  be  proportional  to 
<jh ,  i.e.  hm  =  CFaH,  in  the  high  Froude  number  limit.  Thus,  the  ‘midship  extrapolation  height’  hm  is  taken 


as 


hm  =  min  ( em  2  7r F2 ,  CF  <jh) 


(9a) 
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Figure  2:  Wave  profiles,  predicted  by  Hogner’s  slender-ship  approximation,  along  the  Wigley  hull  at  four 
Froude  numbers  F  =  0.15,  0.25,  0.35  and  0.5.  The  curves  marked  h  =  hm  or  h  =  h[>  correspond  to  an 
extrapolation  height  h  that  is  constant  for  —0.5  <x<  0.5  and  equal  to  the  basic  extrapolation  heights  hm  or 
hb  defined  by  (10a).  The  curves  marked  A  =  1  correspond  to  fcoo  given  by  (11)  and  A  =  1  (no  parabolic  filter) 
in  (7). 

We  then  have  hm  =  em27 rF2  or  hm  =  C^<jh  for  F  <  Fm  or  Fm  <  F  where  the  transition  Froude  number 
Fm  between  the  low-speed  and  high-speed  regimes  is  given  by 


F=0.35 
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J  o  \ 

if 
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F  = 
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Numerical  tests  for  the  Wigley  hull  show  that  the  physically  reasonable  choices  em  ~  0.06  and  ~  1  are 
satisfactory.  These  values  of  em  and  yield  a  transition  Froude  number  Fm  «  0.4  (for  which  we  have 
Ao  ~  1)  for  a  typical  hull  slenderness  aH  =  0.06. 

Ship  bow  waves  and  stern  waves  —  in  the  vicinity  of  a  ship  bow  and  stern  where  the  hull  geometry  varies 
rapidly  —  differ  significantly  from  a  transverse  wave  with  wavelength  approximately  equal  to  Ao.  In  particular, 
a  ship  bow  wave  is  significantly  shorter,  higher,  and  consequently  sharper  than  waves  along  a  ship  hull  aft 
of  the  bow  wave.  A  ship  bow  wave  also  contains  more  short  divergent  waves  and  is  more  directly  affected 
by  the  hull  geometry,  specifically  the  shape  of  the  ship  bow.  Indeed,  a  ship  bow  wave  is  strongly  affected  by 
nearfield  effects  associated  with  a  sudden  change  of  flow  direction  at  a  ship  bow;  [11].  In  particular,  while 
waves  along  a  ship  hull  aft  of  the  ship  bow  are  nearly  sinusoidal,  with  wavelength  approximately  equal  to  Ao  as 
already  noted,  the  portion  of  the  wave  profile  between  a  ship  stem  and  the  crest  of  the  bow  wave,  i.e.  the  bow 
wave  front,  is  well  approximated  by  a  parabola  [11,12]  .  The  base  of  this  parabola,  i.e.  the  (nondimensional) 
distance  lfr0nt  =  Lfront  / Ls  from  the  ship  stem  to  the  a:— location  of  the  bow  wave  crest,  is  shown  in  [11,13] 
to  be  approximately  given  by 


tfront  =  CF2/(1  +  Vs/^D)  »  CF2/(1  +  Fbd )  =  CF2/(  1  FF/Va^) 
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Figure  3:  Sinkage  (top  left  corner),  trim  (top  right),  nearfield  drag  (bottom  left)  and  farfield  drag  (bottom 
right)  predicted  by  Hogner’s  slender-ship  approximation  for  the  Wigley  hull  in  the  range  0.15  <  F  <  0.5.  The 
curves  marked  h  =  hm  or  h  =  hb  correspond  to  an  extrapolation  height  h  that  is  constant  for  —0.5  <  x  <  0.5 
and  equal  to  the  basic  extrapolation  heights  hm  or  h.b  defined  by  (10a).  The  curves  marked  A  =  1  correspond 
to  fcoo  given  by  (11)  and  A  =  1  (no  parabolic  filter)  in  (7). 


with  a  proportionality  factor  C  ~  1.1  for  a  wedged  shaped  bow.  Fbd  is  the  Froude  number  (1)  associated 
with  the  characteristic  transverse  dimension  LBD ,  which  is  used  here  instead  of  the  ship  draft  D.  Thus,  £  front 
is  0(F2)  as  F— >  0  and  O(F)  as  F— >■  oo.  An  effective  wavelength  Xbow  that  approximately  corresponds  to  the 
length  £  front,  of  the  front  of  the  bow  wave  can  be  taken  as  C' if  ront  with  a  proportionality  factor  4  <  C'  <  5. 
The  extrapolation  height  hb  related  to  the  filtering  of  short  waves  at  a  ship  bow  is  taken  as  a  small  fraction 
0  <  Eb  -C  1  of  the  effective  bow  wavelength  Xbow .  Thus,  hb  is  chosen  as  hb  =  CstEbF2/(  1  +  F/V<jh)  with 
C*  =  C'C  ~  5.  The  value  C*  =  5  is  used  here.  This  expression  for  the  extrapolation  height  hb  becomes 
meaningless  at  high  Froude  numbers,  as  already  noted.  Indeed,  as  for  the  ‘midship  extrapolation  height’  hm, 
the  ‘bow-wave  extrapolation  height’  hb  can  be  assumed  to  be  proportional  to  aH ,  i.e.  hb  =  aH ,  in  the  high 
Froude  number  limit.  Thus,  hb  is  taken  as 


hb 


min 


5  ebF2 
+  F/V^ 


(9b) 


The  transition  Froude  number  Fb  between  the  low-speed  and  high-speed  regimes  is  given  by 


Fh  = 


_ 


10  Eb 


^1+  \Jl+  20eb/CB  ^ 


Numerical  tests  for  the  Wigley  hull  show  that  the  physically  reasonable  choices  eb  ~  0.04  and  CB  ~  0.2  are 
satisfactory.  These  values  of  eb  and  CB  yield  a  transition  Froude  number  Fb  ~  0.4  for  a  typical  hull  slenderness 
<jb  =  0.06. 

The  transition  Froude  numbers  Fm  «  0.4  and  Fb  ss  0.4  are  consistent.  These  transition  Froude  numbers 


are  strictly  equal  if 


8tt £m/C^  =  Al  /(1  +  \J\  +  Ab)2  with  Ab  =  20eb/C b 

This  relation  yields  ~  1  for  em  =  0.06,  eb  =  0.04  and  Cb  =0.2. 

Alternatively,  the  midship  and  bow-wave  extrapolation  heights  hm  and  hb  can  be  expressed  as 


hm  =  2nem  min  ( F2 ,  F2 ) 


hb  =  5  eb  min 


F2 

1+F/V IP* 


(10a) 


As  already  noted,  numerical  tests  for  the  Wigley  hull  show  that  Ft  «  0.4,  em  «  0.06  and  eb  «  0.04  are 
satisfactory  choices,  and  expressions  (10a)  are  then  used  here  with  Ft  =  0.4,  em  =  0.06  and  eb  =  0.04. 

The  hull  geometry  near  the  ship  stern  typically  varies  rapidly,  as  in  the  bow  region.  Thus,  we  can  expect 
that  the  extrapolation  height  h  should  decrease  near  the  stern,  as  in  the  bow  region.  However,  the  wave 
generated  by  a  ship  stern  is  not  as  well  understood  as  the  bow  wave  —  for  which  considerable  information 
is  known,  e.g.  [11-14]  —  and  is  more  complicated  due  to  viscous  effects.  For  lack  of  greater  knowledge,  the 
‘stern  extrapolation  height’  hs  is  then  simply  taken  proportional  to  the  midship  extrapolation  height  hm,  i.e. 


h,  =  C™  h 


m 


(10b) 


Numerical  tests  for  the  Wigley  hull  show  that  C ™  =  0.7  is  a  reasonable  choice.  However,  for  a  ship  hull  with 
a  transom  stern  for  which  the  flow  extends  smoothly  aft  of  the  transom,  the  geometry  of  the  ship  hull  surface 
and  related  wake  model  aft  of  the  transom  does  not  vary  rapidly  near  the  stern,  and  the  choice  C™  «  1  is 
appropriate. 

Expressions  (10a)  and  (10b)  yield  hb  <  h.  Thus,  in  (6)  and  the  related  cutoff  wavenumber  k ^  in  (7) 
can  be  taken  as  hoo  =  hb  and  koo  —  \JhF2/hb  —  1  with  hb  given  by  (10a).  We  then  have 


k 


OO 


if  F  <  Ft 


koo 


F  ll  +  Ft/VaV-£bFt2/F 2 
Ft\  £b 


if  Ft  <  F 


(11) 


These  relations  show  that  koo  increases  with  respect  to  the  Froude  number  F,  approximately  in  proportion  to 
F  or  y/F  in  the  high-speed  or  low-speed  regimes  Ft  <  F  and  F  <  Ft,  respectively.  The  smallest  value  of  k^ , 
attained  in  the  limit  F  =  0,  is  y/1  —  £b /y/£b  ~  1  l\f£\>-  We  then  have  5  <  koo  for  £b  =  0.04. 

The  nondimensional  wave  drag  coefficient  Cwave  =  Dwave/(pV2Ll)  that  is  associated  with  the  energy 
transported  by  the  farfield  waves  is  determined  in  terms  of  the  ‘farfield’  wave  spectrum  function  Slfar  = 
S(k;  x  =  x8)  evaluated  at  the  ship  stern  x  =  xs.  Specifically,  the  farfield  wave  drag  coefficient  C^v e  is  defined 
by  the  classic  Havelock  formula  as 

k^ar 

Cwave  =  r-  /  dk  e~5  (fe/fefar)4  y/1  +  k*  (S2  +  Sf )  with  Sr  =  Re(S'far)  =  Im(5far)  (12) 

Z7T  J  —  k  far 


where  the  filter  function  e~5  ^k/k  ^  is  inserted  to  ensure  more  robust  convergence.  The  truncation  wavenum- 
ber  fcfar  is  related  to  the  wave  spectrum  function  at  the  ship  stern  x  =  xa  as  already  noted,  and  can  then  be 
taken  as  fcfar  =  \J§  F2/hs  —  1.  Expressions  (10a)  and  (10b)  yield 


fcfar 


/ 


2.5/7 r 

^  s  t m 


-l 


if  F  <  Ft 


kial 


F  I  2.5/tt 

Ft\j  C™£m 


F? 

F2 


if  Ft  <  F 


(13) 


These  relations  show  that  fcfar  is  independent  of  the  Froude  number  (approximately  equal  to  4.2  for  £m  =  0.06 
and  C™  =  0.7)  in  the  low-speed  regime  F  <  Ft,  and  increases  approximately  in  proportion  to  F  in  the 
high-speed  regime  Ft  <  F.  The  low-speed  relation  (13)  is  consistent  with  the  relation  given  in  [3]  . 

The  variations  of  the  midship  and  bow- wave  extrapolation  heights  hm  and  hb  and  of  the  cutoff  wavenumbers 
koo  and  fcfar  with  respect  to  the  Froude  number  F  within  the  range  0  <  F  <  0.8  are  depicted  on  the  left  and 
right  sides  of  Fig.l,  respectively,  for  Ft  =  0.4,  crH  =  0.06,  £b  =  0.04,  £m  =  0.06  and  C™  =  0.7. 

For  illustration  purposes,  the  Hogner  slender-ship  approximation  is  considered  for  the  classic  Wigley 
parabolic  hull.  Fig. 2  depicts  wave  profiles  along  the  Wigley  hull  at  four  Froude  numbers  F  =  0.15,  0.25, 
0.35  and  0.5.  Three  cases  are  considered  in  Fig. 2.  Two  of  the  cases,  marked  h  =  hm  or  h  =  hb,  correspond  to 
an  extrapolation  height  h  that  is  constant  for  —0.5  <x<  0.5  and  equal  to  the  basic  extrapolation  heights  hm 


Figure  4:  The  left  side  depicts  the  functions  xm  =  nmF2  and  xb  =  fibF2/ (1+Fbd)  with  jirn  =  5,  /j,b  =4  and 
crH  =  0.06.  The  right  side  depicts  the  variation  of  the  extrapolation  height  h(x)  defined  by  (14a)-(14d)  with 
(10a)  and  (10b)  for  F  =  0.15,  0.25,  0.35,  0.5  and  Ft  =  0.4,  aH  =  0.06,  £b  =  0.04,  ft b  =  4,  =  5,  em  =  0.06, 

C™  =  0.7  and  ls  =  0.4. 

or  hb  defined  by  (10a).  The  third  case  in  Fig. 2,  marked  A  =  1,  corresponds  to  koo  given  by  (11)  and  A  =  1 
(parabolic  filter  turned  off)  in  (7).  Thus,  the  cutoff  wavenumber  koo  is  the  same  for  the  curves  marked  h  =  hb 
or  A  =  1 ,  and  differences  between  these  curves  illustrate  both  the  practical  necessity  of  removing  short  waves 
and  the  effectiveness  of  the  parabolic  filter  (8)  and  the  related  filter  function  A  in  the  Fourier  integral  (7). 
The  wave  profiles  obtained  using  h  =  hm  or  h  =  hb  differ  significantly.  In  particular,  the  choice  h  =  hm  yields 
significantly  lower  bow  waves  than  the  choice  h  =  hb,  in  accordance  with  the  property  that  a  ship  bow  wave  is 
shorter,  higher,  and  more  influenced  by  short  waves  and  nearfield  effects  associated  with  the  rapid  variation  of 
the  hull  geometry  at  the  bow,  than  waves  aft  of  a  bow  wave.  Differences  between  the  curves  h  =  hm  or  h  =  hb 
are  significant  also  aft  of  the  bow  wave,  notably  near  the  ship  stern  x  =  —1/2  and  at  high  Froude  numbers. 
In  particular,  the  curves  h  =  hb  contain  oscillations  (most  clearly  apparent  for  F  =  0.25  and  F  =  0.35)  that 
do  not  appear  in  the  curves  h  =  hm.  The  differences  between  the  curves  marked  h  =  hm  or  h  =  hb  illustrate 
that  h  =  hb  or  h  =  hm  are  better  choices  for  a  bow  wave  or  waves  aft  of  a  bow  wave,  respectively,  as  already 
explained. 

Fig. 3  shows  the  sinkage  and  trim  (top  left  and  right  corners),  and  the  nearfield  and  farfield  wave  drags 
(bottom  left  and  right  corners)  —  obtained  via  integration  of  the  pressure  at  the  hull  surface  or  the  Havelock 
formula  (12),  respectively  —  for  0.15  <  F  <  0.5,  that  correspond  to  the  curves  marked  h  =  hm,  h  =  hb 
or  A  =  1  considered  in  Fig. 2.  The  nondimensional  sinkage,  trim,  and  wave  drag  depicted  in  Fig. 3,  and  in 
subsequent  figures,  are  defined  as  in  [3] .  The  extrapolation  heights  h  =  hm  or  h  =  hb  yield  nearly  identical 
predictions  of  trim  and  farfield  drag.  For  the  sinkage  (top  left  corner),  differences  between  the  curves  h  =  hm 
or  h  =  hb  are  small,  and  increase  as  the  Froude  number  F  increases.  Differences  for  the  nearfield  drag  (bottom 
left  corner)  are  also  relatively  small,  although  larger,  and  likewise  increase  as  F  increases.  The  choice  h  =  hm 
yields  a  lower  nearfield  drag  than  the  choice  h  =  hb,  in  accordance  with  the  property  illustrated  in  Fig. 2  that 
h  =  hm  yields  a  lower  bow  wave  than  h  =  hb-  Whereas  Fig. 2  shows  relatively  large  differences  between  the 
wave  profiles  for  h  =  hb  or  A  =  1,  the  curves  h  =  hb  or  A  =  1  in  Fig. 3  can  hardly  be  distinguished.  This 
result  justifies  the  removal  of  short  waves.  Indeed,  short  waves  have  much  smaller  effect  on  integrated  flow 
quantities  (sinkage,  trim  and  wave  drag)  of  greatest  practical  importance  than  on  the  wave  profiles  along  the 
ship  hull  and  the  pressure  distribution  at  the  hull. 


4  Bow-to-stern  variation  of  extrapolation  height 

The  foregoing  analysis  shows  that  the  extrapolation  height  h  in  (8)  should  be  allowed  to  vary  along  the  ship 
hull.  The  variation  of  h  along  the  hull  is  then  taken  as 

h  =  hm  (hm  hb )  Bb  (x)  ( hjyib  hs )  Bs  (x)  with  hmb  =  hm  ( hm  hb )  Bb(x  =  1/2)  (14a) 

The  functions  Bb{x)  and  Bs(x)  are  equal  to  1  at  the  bow  x  =  1/2  or  the  stern  x  =  —1/2,  respectively.  The 
length  of  the  bow  region  where  we  have  Bb  ~  1  is  related  to  the  distance  £  front  ~  F2/{  1  +FBD)  between  the 
ship  stem  and  the  crest  of  the  bow  wave,  which  defines  the  front  of  the  bow  wave.  We  then  impose  Bb  ~  1  for 
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Figure  5:  Wave  profiles,  predicted  by  Hogner’s  slender-ship  approximation,  along  the  Wigley  hull  at  four 
Froude  numbers  F  =  0.15,  0.25,  0.35  and  0.5.  The  curves  marked  h  =  hm  or  h  =  hb  correspond  to  an 
extrapolation  height  h  that  is  constant  for  —0.5  <  x  <  0.5  and  equal  to  the  basic  extrapolation  heights  hm 
or  hb  defined  by  (10a).  The  curves  h  =  hv  correspond  to  an  extrapolation  height  h.(x)  that  varies  from  the 
bow  to  the  stern  as  defined  by  (14a)-(14d)  with  (10a)  and  (10b).  The  vertical  lines  mark  the  boundaries  of 
the  transition  region  defined  by  (14b). 


1/2  —  /ibF2 / (1  +  Fbd)  <  x  <  1/2.  The  parameter  0  <  fib  controls  the  extent  of  the  bow  region  where  Bb  ~  1. 
We  also  impose  Bb  «  0  at  a  distance  aft  of  the  ship  bow  that  is  related  to  the  distance  between  the  ship  stem 
and  the  intersection  of  the  bow  wave  with  the  mean  free  surface.  This  distance  (the  length  of  the  bow  wave) 
is  shown  in  [13]  to  be  approximately  given  by  2.3 F2 .  We  then  impose  Bb  ~  0  for  x  <  1/2  —  fimF2 .  The 
parameter  0  <  fim  controls  the  extent  of  the  region  within  which  the  function  Bb(x)  is  not  approximately  nil. 
Thus,  the  function  Bbix)  varies  from  1  to  0,  and  the  related  extrapolation  height  h  defined  by  (14a)  varies 
from  hb  to  hmi  within  the  region 


l/2~  fimF2<x<  1/2- fibF2/{l  +  FBD) 


(14b) 


The  length  [fim  —  fib /(1+Fbd)]F2  of  this  transition  region  is  controlled  by  the  parameters  fib  and  fim,  and  is 
positive  if  fib  / fim  <  1+Fbd  =  1 +F/V crH .  The  distance  HmF 2  related  to  the  length  of  the  bow  wave  is  smaller 
or  greater  than  1  in  the  low  or  high  Froude  number  regimes  F  <  1/ y/Hm  or  1/ y/iim  <  F ,  for  which  expression 
(14a)  for  hmb  yields  bmb  =  hrn  or  hmb  <  hm,  respectively.  Values  of  ^  and  fim  within  the  ranges  3  <  /it,  <  5 
and  4  <  /im  <  6  are  reasonable.  Numerical  tests  show  that  =  4  and  /rm  =  5  are  satisfactory  choices.  The 
choice  nm  =  5  yields  1  /^/JFn  =  1  for  F  m  0.45.  The  functions  xm  =  /xmF 2  and  xb  =  fib F2/ ( 1  +  Fbd )  are 
depicted  on  the  left  side  of  Fig. 4  for  fim  =5,  fib  =4  and  aH  =  0.06.  We  have 
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Figure  6:  Sinkage  (top  left  corner),  trim  (top  right),  nearfield  drag  (bottom  left),  farfield  drag  (bottom  right) 
predicted  by  Hogner’s  slender-ship  approximation  for  the  Wigley  hulHn  the  range  0.15  <  F  <  0.5.  The 
curves  marked  h  =  hm  or  h  =  hb  correspond  to  an  extrapolation  height  h  that  is  constant  for  —0.5  <  x  <  0.5 
and  equal  to  the  basic  extrapolation  heights  hm  or  hb  defined  by  (10a).  The  curves  h  =  hv  correspond  to 
an  extrapolation  height  h(x)  that  varies  from  the  bow  to  the  stern  as  defined  by  (14a)-(14d)  with  (10a)  and 
(10b). 


A  continuous  transition  between  hb  and  hm  can  be  defined.  For  instance,  the  transition  function 

2  x  —  1  +  F2  [  /rm  +  Hb /(l  +  Fbd)  ] 


2 Bb(x)  =  1  +  tanh  : 


with 


-1  _  1 

—  <  x  <  - 

2  ~  ~  2 


(14c) 


F2[nm  -Vb/{l  +  FBD)]/2 

and  Hb  =  4,  fim  =  5  yields  Bb  ~  0.98  for  x  =  1/2  —  jib  F2/ ( 1  +  Fb d )  and  Bb  «  0.02  for  x  =  1/2  —  nmF2  . 

The  function  Bs  associated  with  the  decrease  of  the  extrapolation  height  h  near  the  ship  stern  is  taken  as 


Bs(x)  =  e-^ix  +  W/F  with  —  1/2  <x<  1/2 


(14d) 


The  parameter  £s  defines  the  stern  region  within  which  h  may  be  decreased  if  the  hull  geometry  changes 
rapidly.  The  value  is  =  0.4  is  a  reasonable  choice  and  is  used  here. 

The  right  side  of  Fig. 4  depicts  the  variation  of  the  extrapolation  height  h{x)  defined  by  (14a)-(14d)  with 
(10a)  and  (10b)  for  F  =  0.15,  0.25,  0.35,  0.5  and  Ft  =  0.4,  aH  =  0.06,  Eb  =  0.04,  ^ ib  =  4,  //m  =  5,  em  =  0.06, 
C™  =  0.7  and  iB  =  0.4.  In  accordance  with  expressions  (14a)-(14d),  Fig. 4  shows  that  the  extrapolation  height 
h(x)  is  smallest  within  the  bow  region,  increases  from  the  bow  value  h  =  hb  to  the  midship  value  h  =  hm 
within  the  transition  region  (14b),  and  decreases  from  the  midship  value  h  =  hm  to  the  stern  value  h  =  hs 
near  the  stern  x  =  —1/2.  Fig. 4  also  illustrates  how  the  values  hb,  hm  and  hs,  and  the  extents  of  the  bow- wave 
and  midship  regions,  and  of  the  transition  region  between  these  two  regions,  vary  with  the  Froude  number  F. 

For  illustration  purposes,  the  Hogner  slender-ship  approximation  is  considered  again  for  the  Wigley  parabolic 
hull  in  Fig. 5.  This  figure  depicts  the  wave  profiles  along  the  Wigley  hull  for  four  Froude  numbers  F  =  0.15, 
0.25,  0.35  and  0.5.  Three  cases  are  considered  in  Fig. 5.  Two  of  the  cases,  marked  h  =  hrn  or  h  =  hb,  corre¬ 
spond  to  an  extrapolation  height  h  that  is  constant  for  —0.5  <  x  <  0.5  and  equal  to  the  basic  extrapolation 
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Figure  7:  Wave  profiles  along  the  Wigley  hull,  at  four  Froude  numbers  F  =  0.15,  0.25,  0.35,  0.5,  predicted  by 
the  NM  theory  with  several  values  of  k l  . 

heights  hm  or  hb  defined  by  (10a).  The  third  case,  marked  h  =  hv,  corresponds  to  an  extrapolation  height  h(x) 
that  varies  from  the  bow  to  the  stern  as  defined  by  (14a)-(14d)  with  (10a)  and  (10b).  The  boundaries  of  the 
transition  region  defined  by  (14b)  are  marked  by  two  vertical  lines  in  Fig. 5.  The  extent  of  the  transition  region 
increases  as  F  increases,  as  already  illustrated  on  the  left  side  of  Fig. 4.  For  F  =  0.15,  the  transition  region 
is  narrow,  whereas  only  one  vertical  line  appears  in  the  bottom  right  corner  of  Fig. 5  because  the  transition 
region  extends  beyond  the  ship  stern  x  =  —0.5  for  F  =  0.5.  In  the  bow  region  to  the  right  of  the  transition 
region,  the  wave  profiles  h  =  hv  and  h  =  hb  cannot  be  distinguished,  and  are  appreciably  higher  than  the 
profiles  h  =  hm  as  already  noted  in  Fig. 2.  To  the  left  of  the  transition  region,  the  profiles  h  =  hv  and  h  =  hm 
are  nearly  identical,  except  near  the  stern  x  =  —0.5  where  the  wave  profiles  h  =  hv  are  lower  than  the  profiles 
h  =  hm  in  accordance  with  (10b).  Within  the  transition  regions  marked  by  vertical  lines  in  Fig. 5,  the  wave 
profiles  h  =  hv  smoothly  merge  between  the  profiles  h.  =  hb  and  the  profiles  h  =  hm  as  expected  from  (14a). 

Fig. 6  shows  the  sinkage  (top  left  corner),  the  trim  (top  right),  and  the  nearfield  (bottom  left)  and  farfield 
(bottom  right)  wave  drags  —  predicted  by  the  Hogner  approximation  for  the  Wigley  hull  -  that  correspond  to 
the  wave  profiles  shown  in  Fig. 5.  In  accordance  with  Fig. 3,  the  curves  h  =  hmi  h  =  hb  and  h  =  hv  can  hardly 
be  distinguished  for  the  trim  and  the  farfield  wave  drag,  and  differences  are  relatively  small  for  the  nearfield 
drag  and,  especially,  the  sinkage.  In  particular,  the  curves  h  =  hv  are  h  =  hb  are  quite  close,  especially  at  low 
Froude  numbers  and  for  the  sinkage. 


5  Neumann-Michell  theory 

The  Neumann-Michell  (NM)  theory  of  ship  waves  given  in  [3]  expresses  the  amplitude  A  of  the  elementary 
wave  E  in  the  wave-spectrum  function  (3)  as  the  Hogner  amplitude  A  =  AH  given  by  (4)  and  a  correction 
Ax .  Specifically,  the  NM  amplitude  function  is  given  by 


A  =  Anm  =  Ah  -f  A^  with  A ^  =  ( kvy  +  i  \J\  +  k2  vz )  (f>t  +  nx  ( k  vz  —  i  \/l  +  k2  vv )  <pd  (15) 
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Figure  8:  Sinkage  (top  left  comer),  trim  (top  right),  nearfield  drag  (bottom  left),  farfield  drag  (bottom  right) 
for  the  Wigley  hull,  in  the  range  0.15  <  F  <  0.5,  predicted  by  the  NM  theory  with  several  values  of  kL- 


Here,  (vv,vz)  =  (riv ,nz)/is  with  is  =  yj ( ny )2  +  (n2)2 .  Furthermore,  <f>d  =  d<j)/dd  and  4>t  =  dcjs/dt  are  the  flow 
velocities  along  the  orthogonal  unit  vectors  d  =  (0,  —  vz,  isv)  and  t  =  (is,  —  nxisv ,  —  nxuz).  The  unit  vectors  d 
and  t  are  tangent  to  the  ship  hull  surface  T,H .  The  unit  vector  d  points  upward  on  the  positive  side  of  Eff, 
where  0  <  y  and  0  <nv ,  or  downward  on  the  negative  side  of  T,H  where  y  <  0  and  ny  <  0.  The  unit  vector  t 
points  toward  the  ship  bow  on  both  sides  of  the  hull  surface  T,H . 

The  Hogner  amplitude  AH  =  nx  is  independent  of  k  and  consequently  does  not  grow  as  k  — >  ±oo.  However, 
the  amplitude  function  A ^  of  the  NM  correction  is  0(k)  as  k  — >  ±oo.  The  NM  wave  potential  associated  with 
the  amplitude  function  A ^  in  the  wave-spectrum  function  S  is  then  influenced  by  short  waves  to  a  greater 
extent  than  the  Hogner  wave  potential  associated  with  the  amplitude  function  AH .  The  short-wave  content 
of  the  NM  wave  correction  potential  can  be  attenuated  by  replacing  k  in  (15)  by  the  function  fc*  defined  as 

A:*  =  k  for  —  kL  <  k  <  kL  k*  =  Aisign(fc)  for  k l  <  |fc|  (16) 

This  modification  does  not  affect  long  waves  associated  with  the  range  —  kL  <  k  <  k^  but  affects  short  waves 
outside  this  long-wave  range. 

Fig. 7  shows  the  wave  profiles  along  the  Wigley  hull,  at  F  =  0.15,  0.25,  0.35  and  0.5,  predicted  by  the  NM 
theory  with  several  values  of  kL-  The  wave  profiles  obtained  for  =  1  and  3  cannot  be  distinguished  for 
F  =  0.15  and  are  fairly  close  for  F  =  0.25.  For  F  =  0.35,  the  wave  profiles  for  kL  =  2  and  3  can  hardly  be 
distinguished,  but  these  profiles  differ  appreciably  from  the  profile  obtained  for  =  1.  For  F  =  0.5,  the  wave 
profiles  are  significantly  influenced  by  kL  -  In  particular,  the  wave  profile  kL  =  1  is  significantly  less  oscillatory 
than  the  profiles  =  2,3  and  5.  The  increasing  effect  of  kL  for  larger  Froude  numbers  illustrated  in  Fig. 7 
is  consistent  with  Fig.l,  which  shows  that  short  waves  become  increasingly  more  important  as  the  Froude 
number  increases.  Fig. 8  shows  that  differences  in  the  sinkage,  trim  and  nearfield  drag  obtained  for  kL  =  1 
and  kL  =  5  are  small.  The  farfield  drag  obtained  for  kL  =  1  is  appreciably  lower,  especially  at  high  Froude 
numbers,  than  the  drags  obtained  for  kL  =  2  or  5,  which  are  relatively  close.  Fig. 8  and  Fig. 7  suggest  that  the 
physically  reasonable  choice  kL  =  2,  which  means  that  waves  with  wavelength  Ao/5  <  A  <  Ao  are  unaffected 
by  the  modification  (16),  is  satisfactory,  and  this  choice  is  used  here. 


Figure  9:  Wave  profiles,  predicted  by  the  NM  theory,  along  the  Wigley  hull  at  four  Froude  numbers  F  =  0.15, 
0.25,  0.35  and  0.5.  The  curves  marked  h  =  hm ,  k^  =  2  correspond  to  an  extrapolation  height  h  that  is 
constant  for  —0.5  <  x  <  0.5  and  equal  to  the  midship  extrapolation  height  hm  given  by  (10a).  The  curves 
marked  h  =  hv ,  k^  =  2  correspond  to  an  extrapolation  height  h(x)  that  varies  from  the  bow  to  the  stern 
as  specified  by  (14a)-(14d)  with  (10a)  and  (10b).  The  value  of  in  (16)  is  taken  equal  to  2  for  the  curves 
marked  h  =  hm,  k^  =  2  or  h  =  hv ,  k^  =  2 .  The  curves  marked  A  =  1,  k^  =  5  correspond  to  k^  =  5,  A  =  1 
(no  parabolic  filter)  and  fcoo  given  by  (11)  in  (7),  i.e.  short  waves  are  not  filtered  in  the  curves  A  =  1,  =  5. 

The  need  for  filtering  short  waves  is  further  illustrated  in  Fig. 9,  which  shows  wave  profiles  predicted  by 
the  NM  theory  along  the  Wigley  hull  at  four  Froude  numbers  F  =  0.15,  0.25,  0.35  and  0.5,  and  in  Fig. 10 
which  depicts  the  sinkage  (top  left  corner),  trim  (top  right),  nearfield  drag  (bottom  left)  and  farfield  drag 
(bottom  right)  predicted  by  the  NM  theory  for  the  Wigley  hull  in  the  range  0.15  <  F  <  0.5.  Three  cases  are 
considered  in  these  two  figures.  The  curves  marked  h  =  hm,  k^  =  2  correspond  to  an  extrapolation  height  h 
that  is  constant  for  —0.5  <  x  <  0.5  and  equal  to  the  midship  extrapolation  height  hm  given  by  (10a).  The 
curves  marked  h  =  hv ,  k^  =  2  correspond  to  an  extrapolation  height  h(x)  that  varies  from  the  bow  to  the 
stern  in  accordance  with  (14a)-(14d)  with  (10a)  and  (10b).  The  value  of  k l  in  (16)  is  taken  equal  to  2  for 
the  curves  marked  h  =  hm ,  &£  =  2  or  h  =  hv ,  kz,  =  2 .  The  main  difference  between  the  curves  h  =  hm  and 
h  =  hv  is  that  the  short  waves  contained  in  the  bow  wave  are  included  in  the  curves  h  =  hv  but  are  ignored  in 
the  curves  h  =  hm .  The  curves  marked  A  =  1 ,  =  5  correspond  to  /cl  =  5 ,  A  =  1  (no  parabolic  filter)  and 

koo  given  by  (11)  in  (7).  Thus,  short  waves  are  not  filtered  in  the  curves  A  =  1,  /cl  =  5.  Fig. 9  shows  that  the 
wave  profiles  A  =  1,  k^  =  5  contain  large  unrealistic  oscillations,  especially  aft  of  the  bow  wave,  and  that  the 
wave  profiles  h  =  hm,  kL  =  2  exhibit  overly  damped  bow  waves,  notably  for  F  =  0.35  and  0.40.  Fig. 10  shows 
that  differences  in  the  sinkage,  trim  and  wave  drag  associated  with  the  alternative  choices  h  =  hm  or  h.  =  hv 
are  appreciable  for  the  wave  drag,  but  insignificant  for  the  sinkage  and  trim.  The  ‘no-filtering’  choice  A  =  1 
and  fcz,  =  5  yields  overly  oscillatory  wave-drag  curves  that  are  significantly  different  from  the  curves  h  =  hm 
or  h  =  hv . 
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Figure  10:  Sinkage  (top  left  corner),  trim  (top  right),  nearfield  drag  (bottom  left),  farfield  drag  (bottom  right) 
predicted  by  the  NM  theory  for  the  Wigley  hull  in  the  range  0.15  <  F  <  0.5.  The  curves  marked  h  =  hm, 
Ul  =  2  correspond  to  an  extrapolation  height  h  that  is  constant  for  —0.5  <x<  0.5  and  equal  to  the  midship 
extrapolation  height  hm  given  by  (10a).  The  curves  marked  h  =  hv ,  ki,  =  2  correspond  to  an  extrapolation 
height  h(x)  that  varies  from  the  bow  to  the  stern  as  specified  by  (14a)-(14d)  with  (10a)  and  (10b).  The  value 
of  fci  in  (16)  is  taken  equal  to  2  for  the  curves  marked  h  =  hm ,  kh  =  2  or  h  =  hv ,  ki,  =  2 .  The  curves  marked 
A  =  1,  fej,  =  5  correspond  to  k^  =  5,  A  =  1  (no  parabolic  filter)  and  koo  given  by  (11)  in  (7),  i.e.  short  waves 
are  not  filtered  in  the  curves  A  =  1,  k^  =  5. 

6  Conclusion 

A  reasonable  and  practical  approach,  based  on  parabolic  extrapolation  within  a  layer  bordering  the  mean 
free-surface  plane  z  =  0,  for  evaluating  waves  in  the  vicinity  of  the  plane  z  =  0  within  the  framework  of 
the  ‘free-surface  Green  function  potential  flow  theory’  (based  on  a  Green  function  that  satisfies  the  radiation 
condition  and  the  Kelvin-Michell  linearized  boundary  condition  at  the  free  surface)  has  been  given.  This 
parabolic-extrapolation  approach  effectively  filters  short  waves,  as  illustrated  in  Fig. 2  and  Fig. 9.  The  reported 
analysis  and  illustrative  calculations  confirm  the  (well-understood)  need  for  removing  ship  waves  that  are  short 
in  comparison  to  the  ship  length,  or  that  may  be  so  short  as  to  be  affected  by  surface  tension  or  viscosity,  and 
consequently  are  not  relevant  in  practice  or  indeed  physically  unrealistic. 

The  short-wave  filter  based  on  parabolic  extrapolation  amounts  to  modifying  the  basic  Fourier  represen¬ 
tation  (2)  of  ship  waves  as  in  (7).  The  filter  function  A  in  the  modified  Fourier  representation  (7)  is  given 
by  (8)  for  the  parabolic  extrapolation  considered  here.  The  height  h  of  the  extrapolation  layer  associated 
with  the  filter  function  A  depends  on  the  Froude  number  and  the  slenderness  of  the  ship  hull,  and  varies 
from  the  ship  bow  to  the  ship  stern,  as  illustrated  in  Fig. 4.  This  bow-to-stern  variation  accounts  for  the 
fact  that  waves  along  the  ship  hull  aft  of  the  bow  wave  differ  from  the  bow  wave.  Indeed,  the  bow  wave  is 
significantly  higher,  shorter  and  therefore  steeper  than  waves  aft  of  the  bow  wave,  is  significantly  affected  by 
nearfield  effects  related  to  the  rapid  variation  of  the  hull  geometry  at  a  ship  bow,  and  consequently  contains 
more  short  wave  components.  The  extrapolation  height  h  is  determined  explicitly,  by  means  of  expressions 
(14)  and  (10),  in  terms  of  the  Froude  number  F,  the  slenderness  aH  of  the  ship  hull,  and  the  location  x 


along  the  ship  waterline.  These  expressions  involve  six  parameters  that  have  a  clear  physical  meaning.  The 
Neumann-Michell  correction  to  the  Hogner  slender-ship  approximation  involves  an  additional  parameter  that 
restricts  the  influence  of  divergent  waves  and  also  has  a  physical  significance.  Numerical  predictions,  notably 
of  the  sinkage,  trim  and  wave  drag,  are  not  overly  sensitive  to  the  values  chosen  for  these  seven  parameters 
associated  with  the  short-wave  filter  considered  here. 

Comparisons  between  experimental  measurements  and  numerical  predictions  given  by  the  NM  theory  or 
the  related  Hogner  slender-ship  approximation  with  the  short-wave  filter  defined  by  (8),  (10),  (14),  (16)  and 

eb  =  0.04  /I,,  =  4  nm  =  5  em  =  0.06  C"1  =  0.7  ts  =  0.4  kL  =  2 

are  reported  in  [15]  for  the  classic  Wigley  parabolic  hull  and  Series  60  model  . 
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